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1. Introduction to Chiral Ward Identities 

By decoupling the heavy particles such as W bosons, Z bosons and heavy quarks from the 
standard model of electroweak theory, we obtain the low energy effective Hamiltonian relevant to 
particular processes classified as the AS = 1 interaction of primary interest in this paper. 

r 10 

^A5=i = *VrtVZY t \z i (ji) + Ty i (ji)]Q i ( f i) (1.1) 

Here, Gp is the Fermi coupling constant and Vn are elements of the CKM matrix. Note that 
1 = —X t jX u where we define X, = VyV ( *. The z,-(ju) and yi(fl) are the Wilson coefficients at the 
scale pL. The 2,-(ju) are the four fermion operators made of u, d, s quark fields. The Wilson 
coefficients for the current-current operators 2i and 22 are of £?(1), those for the QCD penguin 
operators (£?;, i = 3,4,5,6) are of &{tt s ), whereas those for the electroweak penguin operators (2,-, 
i = 7,8,9, 10) are of 0(a). 

In this paper, we focus on the electroweak penguin operators Q7 and 2s- 

3 

Qi = ^{s a d a )v-A Yj e q(yp<lp)v+A (1-2) 

~ q=u,d,s 

3 

2s = -(* a rf J 8)y-A Yj e q(qp1u)v+A (1-3) 

According to the group theory analysis, 27 and 2s belong to the (8,8) irreducible representation 
(irrep) of SU (3)t (8) >S£/ (3)r. We can further decompose 2(8,8) i nto A/ = 1/2 and A/ = 3/2 irrep. 

2(8,8) = -(^)v-a(2mm- dd- ss) v +a (1-4) 

2(8.8) 3/2 = (^)v-A(KM)v+A + (jM)v-A(^)7+A-(*rf)v-A(^)v+A (1-5) 
2(8.8) 1/2 = (^)v-a(Sm)v+A- («<)v-a(k^)7+A- (jii)v-A(^)v+A (1-6) 

We can rewrite 2(8.8) as follows: 

n _ ^^=3/2 ^=1/2-1 n - 

2(8,8) - 2^(8,8) +2(8,8) J <- L7 ) 

Using chiral perturbation theory, we can analyze the chiral behavior of hadronic matrix el- 
ements of the 2(8.8) operators. The details of this analysis will be presented in Ref. At the 
lowest order in chiral perturbation theory, the ratio of the K + — ► n + amplitudes satisfies the follow- 
ing chiral Ward identity |g [| [j}. 

W ^ ^575 > limW = 2 (1.8) 

In this paper, one of the main goals is to numerically check the above Ward identity, which serves 
as an important probe to test the validity of using improved staggered fermions in the calculation 
of weak matrix elements. 
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parameter 


value 


gauge action 


Wilson Plaquette 


fermion action 


HYP Staggered 


P 


6.0 (quenched QCD) 


# of confs 


218 


lattice 


16 3 x 64 


quark mass 


0.01,0.02, 0.03, 0.04 



Table 1: Parameters for the numerical study 



2. Numerical Study on Ward identities 

We use an ensemble of 218 gauge configurations whose details are summarized in Tab. |]. 
We use the p meson to set the scale I /a = 1.95 GeV. We also use the kaon to set the physical 
strange quark mass. Details of this procedure are presented in Ref. Throughout this paper, we 
consider only particles composed of degenerate quarks (m s = m£). One of the key ingredients of 
this numerical study is that we use improved staggered fermions made of HYP fat links when we 
construct the operators and calculate the matrix elements. 

2.4 i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 2.2 
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Figure 1: W for Q 7 (left) and Q & (right) 



We measure W of Eq. ( |1.8j ) for four different quark masses. The results for Q-j and Q% are 
presented in Fig. [j]. In Fig. |] we plot the ratio W as a function of m 2 K . We fit the data to the 
following form suggested by quenched chiral perturbation theory []3|, [[]] : 



/(<) 



1 + 3 zr-§TT 1 °g 



(4n ff & V(47T/)2 



+ ci(m K )+c 2 (m K 



2 \2 



(2.1) 



Here, even though we know that Co = 2 from chiral perturbation theory, we fit the data by con- 
sidering cq a free parameter to see how well it agrees with the prediction of chiral perturbation 
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Qi 




CO 


2.1206 ±0.0820 


2.0574 ± 0.0734 


c\ 


4.1970 ±0,2298 


4.6681 ± 0.2194 


C2 


-3.9582 ± 0.4876 


-3.9037 ± 0.4538 



Table 2: Fit parameters of the Ward identity for Qi and Q%. 

theory. The fitting results are summarized in Tab. ||. Note that the cq values are consistent with the 
theoretical prediction of Co = 2 within a statistical uncertainty of 1.5 a for Qj and 0.8 a for 

Taking into account the fact that we neglect finite volume effects, higher order corrections in 
chiral perturbation theory and two-loop corrections in the matching, the agreement is quite good. 
This suggests that we should be able to obtain a solid result for the Qi and Qg matrix elements in 
the chiral limit using improved staggered fermions. 



3. Numerical Study on Individual Amplitudes 



In Ref. [|3j], Laiho and Soni presented the results of the chiral log corrections to the individual 
matrix elements of the 2(8,8) operator using (partially) quenched chiral perturbation theory. For the 
AI = 3/2 amplitude, the chiral behavior is 



(* + iC)V + ) 
f 2 

J 71 



1-2- 



(47T/)- 



■log 



mi 



For the AI = 1/2 amplitude, the chiral behavior is 



sub 



1 + - 



nn 



■log 



{Ajlff 



mi 



{Anff 



+ bi(m K )+b 2 (m K 



2 \2 



+ d\{m K ) +d 2 {m K 



2 \2 



(3.1) 



(3.2) 



Pn ~T'(4*/) 2 
where the subscript sub means a kind of subtraction method, which will be explained later. Here, 
note that the £>, and d[ parameters cannot be determined from the chiral perturbation theory except 
for the ratio do/bo = 2. 

In the case of the AI = 3/2 amplitudes, there is no mixing with unphysical lower dimension op- 
erators. However, this is not the case for the AI = 1/2 amplitudes. There is mixing with unphysical 
lower dimension operators in the next to leading order 0(m\). The mixing coefficients are power 
divergent of 0(1 /a 2 ). Hence, even though they are suppressed by m\, the size of this mixing dom- 
inates the physical signal due to the divergent coefficient in the range of quark masses that we use 
in this numerical study. Therefore it is essential that the mixing with unphysical lower dimension 
operators should be subtracted non-perturbatively. There are many ways to handle this. The de- 
tails for different options will be discussed in Ref. Here, we focus on the best method that we 
choose for this numerical study. Regardless of the alternative methods to calculate (0|<2|g 3 |^°), if 
we take the derivative of the K° —>0 amplitude with respect to the strange quark mass and take the 
limit of m s — > mj, we can select only the divergent mixing coefficient up to the next leading order 
of the quenched chiral perturbation theory. In other words, 

d 



lim 

m s ^m d dm. 



<0\Q 



(1/2), o N 
(8,8) l A / 



C\B Q . 



(3.3) 
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where Bo is defined as 



m 2 K =Bo(m s + m d ) 



(3.4) 



Here, c' 4 is the divergent mixing coefficient in the notation of Ref. [J3]] . The details on the derivation 



of Eq. (3.3) will be explained in Ref. [{[]]. At any rate, we use Eq. (3.3) to determine c\. Using 



this c r 4 , we can subtract the divergent mixing away. Therefore, the final results of the matrix el- 
ement (7l + \Q^^ \K + ) su \, do not contain any contribution from the mixing with unphysical lower 
dimension operators. 
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Figure 2: (n + \Qf & % ) \K + ) / fl matrix elements for g 8 (left: AI = 3/2) and (right: AI = 1 /2) 

We measure the K + — > n + amplitudes of the Q& operator for four different quark masses. 
The results are presented in Fig. |2[ Here, we normalize the amplitudes by /J for the analysis 
convenience without loss of generality. 1 We fit the data to the functional form suggested by the 
quenched chiral perturbation theory which is given in Eq. and Eq. Q. As you can see m 
Fig. ||, the intercepts of the AI = 1 /2 and AI = 3/2 amplitudes respect the Ward identity (do /bo = 2) 
within statistical uncertainty. 



4. Gold-plated Ratio 

In Ref. [Bt], Becirevic and Villadoro proposed that it is possible to classify ratios of the K + — > 
7t + amplitudes as gold-plated and silver-plated. Here, the gold-plated ratio means that it does not 
contain any chiral log correction at the one loop level and correspondingly the leading finite volume 
effects also cancel off. This gold-plated ratio (R) is defined as 



R 



(7t + \Q<i /2) \K + ) 

(K+\Q { J /2) \K+) 



(4.1) 



Note that by itself cannot possess chial logs in the limit of m s = in quenched QCD. 
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Figure 3: R vs. M\ 

We measure R over four quark masses and the results are presented in Fig. ||. The quadratic 
fitting results in Fig. || are 

R{m 2 K ) = 3.255(24) + 0.235(154) (m|) - 1.153(247)(m|) 2 (4.2) 

Here, note that the coefficient of the quadratic term is reasonable (of £?(1)) whereas the linear term 
has a relatively tiny coefficient. This weak dependence on the linear term might be interpreted as 
a signal for the gold-plated ratio. However, the noticeable dependence on the quadratic term gives 
us some doubt on whether we call R gold-plated literally. 
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